Numerical simulation of stochastic motion of vortex loops under action 
of random force. Evidence of the thermodynamic equilibrium. 
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f"^ ■ Numerical simulation of stochastic dynamics of vortex filaments under action of random 

(Langevin) force is fulfilled. Calculations are performed on base of the full Biot-Savart law for 
£SJ , different intensities of the Langevin force. A new algorithm, which is based on consideration of 

crossing lines, is used for vortex reconnection procedure. After some transient period the vortex 
tangle develops into the stationary state characterizing by the developed fluctuations of various 
physical quantities, such as total length, energy etc. We tested this state to learn whether or not it 
the thermodynamic equilibrium is reached. With the use of a special treatment, so called method of 
weighted histograms, we process the distribution energy of the vortex system. The results obtained 
demonstrate that the thermodynamical equilibrium state with the temperature obtained from the 
' (— i ' fluctuation dissipation theorem is really reached. 

q ■ PACS-numbers: 67.40.Vs 98.80.Cq 7.37.+q 

a 

Quantized vortices appeared in quantum fluids and other systems play a fundamental role in the properties of the 
latter. For that reason they have been an object of intensive study for many years (for review and bibliography see e.g. 
[H). The greatest success in investigations of dynamics of quantized vortices has been achieved in relatively simple 
cases such as a vortex array in rotating helium or vortex rings. However these simple cases are rather exception than 
a rule. Due to extremely involved dynamics initially straight lines or rings evolve to form highly entangled chaotic 
structure. Thus, the necessity of statistic methods to describe chaotic vortex loop configurations arises. A most 
tempting way is to treat vortices as a kind of excitations and to use thermodynamic methods. One of first examples 
of that way was an use the Landau criterium for critical velocity where vortex energy and momentum were applied 
to relation having pure thermodynamic sense. More extended examples would be the famous Kosterlitz-Thouless 
theory or its 3D variant intensively being developed currently {for review and bibliography see e.g. 0). In the 
examples above and in many other it is assumed that chaotic vortex configuration is in thermal equilibrium and their 
statistics obeys the Gibbs distribution. That belief is based on fundamental physical principles and can be justified 
in a standard way considering vortex loops as a subsystem submerged into thermostat and exchanging with energy 
with the latter. However numerous experiments on counterflowing Hell, and direct numerical simulations of vortex 
line dynamics 3, 0] convincingly demonstrate that this dynamics is essentially nonequilibrium and possesses all 
features inherent in turbulent phenomena. Thus a question arises how a thermal equilibrium is destroyed and what 
mechanisms are responsible for that. To answer that question we have firstly to understand in details how a thermal 
equilibrium in vortex loop configuration space is established. A general principle of maximum entropy does not give 
any details the dynamic details are absorbed by a temperature definition. It is well known however that the Gibbs 
distribution can be alternatively obtained on the basis of some reduced model like kinetic equations or Fokker-Planck 
equation (FPE). That way of course is not of such great generality as a principle of maximum entropy, but instead it 
allows to clarify the mechanisms how the Gibbs distribution established [5| . 

In the presented paper we report preliminary results of numerical study on dynamics of vortex tangle under action 
of random (Langevin) forcing delta correlated both in space and time. The data obtained were tested to learn whether 
or not the stationary state reached in numerical experiment is the thermodynamical equilibrium state. With use a 
special treatment, so called method of weighted histograms, to process distribution energy of the vortex system. The 
results obtained demonstrated that the thermodynamical equilibrium with the temperature obtained from the the 
fluctuation dissipation theorem is really reached. 
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II. THE NUMERICAL SIMULATION AND RESULTS 



We consider the dynamics of vortex loops in three-dimensional space with no boundaries. The equation of motion 
of the vortex line elements is supposed to be: 

s = s B + a(s' x (v n -s B ))-o;V x [s' x (v n -s s )] + f(£, t), (1) 

where is the propagation velocity of the vortex filament at a point s, defined by Biot-Savart low; s(£,t) is the 
radius-vector of the vortex line points; v„ is the normal velocity of the superfluid helium; £ is a label parameter, 
in this case it the arc length; s' is the derivative wrt the arc length, a, a' are the friction coefficients, describing 
interaction of vortex filament with normal component, and f (£, t) is the Langevin force. Further we will take v„ = 
and neglect the term with a! . The Langevin force is supposed to be a white noise with the following correlator 

(fi(ei,*i)f J -(6,t 2 )> = DSyS (ti -t 2 )S (&-&). (2) 

Here i,j are the spatial components; t±, t 2 are the arbitrary time moments; £i, £2 define any points on the vortex line; 
D is the intensity of the Langevin's force. Let us consider the probability distribution functional 5] defined as 

V(M0},t) = {5(s(0-s(U)))- (3) 

The Fokker-Planck equation for the time evolution of quantity V({s(^)},t) can be derived from equation of motion 
(P) in standard way (see e.g. 0,0|) 
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(f(f)f (0> *(£-0 S(h - h) 5 Vum = 0. 

As it was shown in [5| the Fokker-Planck equation (|4|) has a stationary solution 

P({s(0}, *) = exp(-H {s} /k B T). (5) 
Here H {s} is the hamiltonian - the energy of the vortex system expressed via the whole line configuration 

H{ s} = E= /^i d 3 r = ^! f hpd^A (6) 



8tt J J |si - s 2 | 


The " temperature" T enterring equation (jJJ is determined with the help of the fluctuation dissipation theorem 

<Ua,ti) M6,t 2 )> = p ^ m) S(h - t 2 ) S m , m . (7) 

In discrete variant of parametrization of the curve, used in numerical simulation, the delta function <5(£i— £2) is 
changed with 1/a, where a is the step along the curve. Thus intensity D of pumping external Langevin force is 
connected to the " temperature" via relation 



D = 2kTa/{p s n). (8) 

Thus, we demonstated that the set of filament agitating by the random white noise is driven into thermodynamical 
equilibrium with the temperature relating to intensity of the random forcing. The main purpose of the the present 
work is to demonstrate it in the direct numerical simulation. 

Details of numerical simulation were described in our early publication Q. The new algorithm for vortex recon- 
nection processes basing on the consideration of crossing lines is used . We run the calculations with a = 0.098, 
D\ = 9.8 • 10~ 5 cm 2 /s, D 2 = 10.2 ■ 10 -5 cm 2 /s. In our calculations we start with an initial vortex configuration 
of twenty four vortex rings (see Fig. 1 (a)). The initial condition was chosen to make the total momentum of the 
system is equal to zero. To check the behavior of the system we had been monitoring the total length of the vortex 
tangle. These quantities for two different intensities of the random forcing were plotted as functions of time in Fig. 
2. One can see, at the beginning the total length rapidly increases. When the vortex tangle becomes dense enough 



FIG. 1: Development of a vortex tangle: t = (o) , t = 6 (6), t = 120 (c) ms. 
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FIG. 2: The lengths of lines as a function of time for different intensity of the Langevin's force D. 



many small vortex loops appear due to reconnection processes (see Fig. 1(b)). These small loops are radiated in the 
ambient space and the total length decrease. Finally a stationary state with the strongly fluctuating is achieved see 
Fig.l (c). We aim now to study this steady state and give some proofs it be in thermal equilibrium. We use for it the 
the weighted histogram analysis method, widely used in the Monte Carlo computer simulations of the Ising model [|[. 

To apply the weighted histogram analysis method we take for every run N configurations of the vortex tangle at 
N different moments in time. We supposed these configurations to be statistically independent. We calculate further 
(With the use of relation ©) the energies of the every configurations. Dividing then the whole interval of energies 
in pieces of width AE we build up the histograms showing relative frequency of meeting the configuration with the 
energies lying in the interval between Ej and Ej + AE (see Fig. 3). Let us introduce the probability density P(Ej) 
of the observing the vortex system state with the energies lying in the interval between Ej and Ej + AE. Obviously, 
the P{Ej) can be simply counted from histogams with the following relation 



Here N(Ej) is the number of states in the interval between Ej and Ej + AE, N = Y^N(Ej) is the full number of 

3 

configurations (for each of two different intensities of the Langevin's force D ) . Consequently assuming that the vortex 
loops system is in the thermal equilibrium we propose that the histograms depicted in Fig. 3 can be describes as well 
with the use of the Gibbs distribution: 



P(Ej) = N(Ej)/(NAE). 



(9) 



P{E 3 ) 



g{Ej)e- E i' kBT 



(10) 



where g(Ej) is the density of states, and T is the temperature (either T\ or T 2 for each of the runs), calculated from 
the fluctuation dissipation theorem O,© for different intensities of the Langevin's force D\, D% correspondingly. 
Relation (JTUJ) included the density of states g(Ej), which is not only known but even is not defined well for set of 
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FIG. 3: Probability density distributions for different intensities of the Langevin's force D which are approximated by Gaussian 
curves. 



continious curves. This problem however can be eliminated with the trick offered in paper [8[. Indeed, it is possible 
to prove that the probability density distribution at Ti can be expressed in terms of the distribution at T\ in the 
following way: 

N ( Ej Ej ) 

Pt (Ei)e (k B T l k BT 2 > 

Pt 2 (E 3 )= Tl[3)e -w. g—. (11) 

AP^P^E^e bTi " bT2 

3 

Using relation (fTTj) we calculated the probability density PT 2 (Ej) for the temperature T 2 via the probability density 
Pti (Ej ) and the compared the result obtained with the initial histogram for the Pr 2 (Ej ) . At this point there appeared 
one difficulty. We mentioned that the discrete variant of the fluctuation dissipation theorem we have to change the 
delta function 8(£±— £2) with quantity 1/a, where a is the step along the curve. But during evolution the distance 
between points does not preserves, the elements of line either shrink or stretch. We nonthclcss retain the initial 
value of the space step in the definition of the temperature but introduce the fitting parameter for it. As a fitting 
parameter we taken half as much again the maximum step along the vortex line using in calculation of dynamics of 
vortex loops. The curves calculated according to equations © (for the temperature T 2 ) and (fTTj) are shown in Fig. 
4. As it can seen the curves are very close to each other. It means that the ensemble of vortex filaments is driven into 
thcrmodynamical equilibrium state. 
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FIG. 4: Probability density distributions at the temperature T2: 1 - the curve calculated according to equation (7); 2 - the 
curve calculated according to equation (8). 

As well it is known that in a thermal equilibrium the variance in the energy (or 'energy fluctuation') is connected 
with the temperature in the following way : 



SE 2 = Ck B T 2 , 



(12) 



■5 



were 

5E^ = ((E-(E)f), (13) 

C = ((E2) — (Ei))/(T2 — T\) is the heat capacity. We determined energy fluctuations according to equation (12): 
8E{T{) = 2.55 • 1CT 16 J, SE(T 2 ) = 2.65 • 10~ 16 J, as well as according to equation (13): 5E{T{) = 2.30 • 10~ 16 J, 
SE(T 2 ) = 2.87 • 1CP 16 J. Obtained energy fluctuations are close to each other. It corroborate that the ensemble of 
vortex filaments is driven into thermodynamical equilibrium state too. 



III. CONCLUSION 

Grounding on results of theoretical study made by one of the authors we propose that the vortex tangle filaments 
undergoing the random "white noise" forcing is driven into thermodynamical equilibrium state. The temperature 
of the vortex system determined by the fluctuation dissipation theorem. With the direct numerical simulation we 
present the proofs of our supposition. The proof was based on the observation that the distribution of the energy 
satisfied to the Gibbs law. It was shown with the method of weighted histograms, widely used approach in statistical 
physics. 

ACKNOWLEDGMENTS 

Authors are grateful to S. Chekmarev for numerous discussions and consultation. This work was partially supported 
by grants 05-08-01375 and 07-02-01124 from the RFBR and grant of the Russian Federation President on the state 
support of leading scientific schools NSH-4366.2008.8. 



[1] Donnelly, R.J. Quantized Vortices in Helium II, Cambridge University Press, 1991. 

[2] G.A.Williams, Vortex-Loop Phase Transitions in Liquid Helium, Cosmic Strings, and High-Tc Superconductors, 

Phys.Rev.Lett., 1999, 82, N6, 1201. 
[3] S. K. Nemirovskii and W. Fiszdon, Chaotic quantized vortices and hydrodynamic processess superfluid helium, Rev. Mod. 

Phys., 1995, 67, Nl, 37. 
[4] A. Chorin, Voticity and Turbulence, Springer- Verlag, New-Yourk, 1994. 

[5] S.K. Nemirovskii, Thermodynamic equilibrium in the system of chaotic quantized vortices in a weakly imperfect Bose gas, 

Teoretical and Matematical Physics, 2004, 141, N 1, 141. 
[6] Jean Zinn- Justin, Quantum Field Theory and Critical Phenomena, Claberson Press, Oxford, 1992. 

[7] LP. Kondaurova and S.K. Nemirovskii, Full Biot-Savart numerical simulation of vortices in He II, J. Low Temperature 
Physics, 2005, 138, N3/4, 555. 

[8] Alan M. Ferrenberg and Robert H. Swendsen, New Monte Carlo technique for studying phase transitions, Phys. Rev. Lett., 
1988, 61, 2635. 



